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'  '  ABSTRACT 

In  the  introduction  of  this  paper  the  basic  principles  of  acousto-optical 
diffraction  in  an  isotropic  medium  are  briefly  reviewed.  Focus  is  on  the 
derivation  of  the  Raman-Nath  equations  for  the  amplitudes  of  the  diffracted  light 
waves  and  on  the  physical  meaning  of  the  various  parameters  occurring  in  this 
diffraction  problem. 

t* 

In  the  first  paragraph,  three  distinct  methods  for  the  numerical  integration 
of  the  truncated  Raman-Nath  system  are  outlined:  Raman-Nath's  elementary  theory, 
Mertens'  perturbation  method,  and  the  N-th  order  approximation  method.  For  each  of 
these  methods,  the  theoretical  results  are  compared  with  experimental  data. 

In  the  second  paragraph^  an  eigenvalue  method  and  an  operational  method  (due 
to  Heaviside-Jeffreys)  are  used  to  integrate  the  truncated  Raman-Nath  system  in 
the  case  of  normal  incidence  of  the  light.  Both  methods  lead  to  closed  form 
expressions  for  the  intensities  of  the  diffracted  light  beams,  which  are  easily 
implemented  on  a  computer. 

A  conqparison  of  the  various  approximation  methods  is  presented,  in  paragraphs 
three  and  four.  ^Fifteen  figures  illustrate  the  validity  of  each  method  in 
comparison  with  experimental  data  obtained  for  different  values  of  the  parameters. 

„  ."''In  paragraph  five,  by  means  of  the  Nth-order  approximation  method,  some  new 
results  are  obtained,  indicating  the  relevance  of  the  three  regime  parameters  in 
distinguishing  between  Raman-Nath  and  Bragg  diffraction. 

In  two  appendices,  examples  are  given  of  the  exact  and  approximate 
integration  of  infinite  and  truncated  Raman-Nath  systems.  ^ 

AMS  (MOS)  Subject  Classifications:  34K99,  78A45;  Physics  Abstracts:  42.10, 

42.80F,K,  43.35,  78.2CH 

Key  Words:  acousto-optics,  diffraction  gratings,  Raman-Nath  systems 
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INTRODUCTION 


Review  of  acousto-optic  diffraction 


The  vibration  of  a  piezo-electric  transducer,  driven  by  a  radio  frequency 
electric  signal,  generates  a  progressive  ultrasonic  wave  into  a  cell  of 
acousto-optic  (AO)  material.  As  sound  passes  through  the  medium,  it  is 
characterized  by  regularly  alternating  compressions  and  rarif actions.  These 
alterations  of  the  internal  strain  cause  a  periodic  variation  of  the  optical 
index  of  refraction  of  the  transparant  medium.  This  type  of  coupling  is 
commonly  called  the  photoelastic  or  elastooptic  effect. 

To  light  impinging  on  the  cell,  the  index  variation  appears  as  a  slowly 
moving  volumetric  phase  grating,  with  spacing  equal  to  the  acoustic  wavelength 
and  grating  depth  proportional  to  the  amplitude  of  the  acoustic  wave. 

Considering  Fig.  1,  where  the  angle  ^  between  the  incident  monochromatic  light 
and  the  acoustic  wavefronts  is  close  to  zero,  we  observe  that  the  emerging 
light  is  diffracted  into  a  set  of  multiple  beams  (orders  at  the  output). 

The  diffracted  beams  are  shifted  up  and  down  in  frequency  by  an  amount  equal 
to  a  multiple  of  the  acoustic  frequency.  The  deflection  angles  of  the  various 
emerging  light  rays  also  depend  on  the  ultrasonic  frequency,  while  their 
intensities  (hence,  the  number  of  output  orders)  are  primarily  determined 
by  the  width  L  of  the  cell  and  the  strength  of  the  sound  field.  The  diffracted 
light  beams  can  be  imaged  on  a  photographic  plate  (using  appropriate  optical 
instrumentation)  so  forming  a  diffraction  spectrum  (Fig.  2). 

To  be  more  precise,  one  has  to  distinguish  between  two  cases,  as  illustrated 
in  Figs.  3a  and  3b.  If  the  input  light  is  incident  perpendicularly  to  a 
rather  narrow  acoustic  cell,  the  diffraction  pattern  shows  a  central  order 
of  undiffracted  light  and  many  side  order  of  diffracted  light.  This  case  of 
multi-order  diffraction  (illustrated  in  Fig.  3a)  is  commonly  called  Raman-Nath 
diffraction  (RN  diffraction).  It  is  achieved  for  rather  intense  sound  fields 
in  the  lower  frequency  range  (<400MHz). 

For  thick  sound  fields  of  higher  frequencies,  essentially  only  two  diffraction 
orders  (zeroth  and  first)  will  exist,  provided  the  input  light  intercepts 
the  sound  cell  at  the  so-called  Bragg  angle,  as  shown  in  Fig.  3b.  This  case 
of  two-order  diffraction  (referred  to  as  Bragg  diffraction)  is  applied  in 
a  majority  of  commercial  acousto-optics  devices  (Bragg  cells). 

The  fact  that  acousto-optic  devices  (either  operating  in  the  Raman-Nath  or 
in  the  Bragg  diffraction  regime)  have  the  capability  of  altering  the  position, 
frequency,  amplitude  and  phase  of  a  (laser)  light  beam,  has  resulted  in  a  variety 


of  applications  in  signal  processing  and  optical  communication. 

Before  proceeding  to  a  brief  discussion  of  the  relative  amplitudes  of  the 
diffracted  light  beams,  let  us  examine  their  directions  of  propagation  and 
frequencies  more  closely. 

The  easiest  way  to  derive  the  deflection  angles  and  the  frequencies  of  the 
scattered  light  beams,  is  to  interpret  the  electromagnetic-acoustic  inter¬ 
action  as  a  photon-phonon  collision  process  [8,13,23].  Recalling  two 
elementary  concepts  of  quantum  wave  mechanics,  i.e.  frequency  is  proportional 
to  energy,  and  wave  number  is  proportional  to  momentum,  we  express  conservation 
of  energy  and  momentum  for  the  interaction  of  the  particles.  Generally  taking 
into  account  multiple  interactions,  reveals  that  the  scattered  light  must 
emerge  in  a  spectrum  whose  frequencies  and  wavenumbers  are  determined  by  the 
conditions  : 


to  =  uj  -  nu8  , 
n 


k  -  1c  -  nlc"  , 
n 


n  =  0,  ±1 ,  ±2 . 


Here  u  and  k  refer  to  the  incident  light  wave:  to  and  k  to  the  n-th  order 

6  n  n 

emerging  light  wave,  while  the  sound  has  an  angular  frequency  gja  and  a 
wavevector  k*.  Note  that  all  symbols  refer  to  the  characteristics  inside 
the  acoustic  medium. 

Restricting  ourselves  to  the  case  of  one  diffracted  order  (n=l),  and  assuming 


exact  momentum  matching  and  downshifted  interaction,  the  corresponding  wave 
vector  diagram  is  easily  drawn  (Fig.  4).  Since  (k*|  <<  (k|,  this  triangle 


is  isosceles,  i.e.  |k^|s 
the  relation 


|  k  | .  It  follows  that  8^ 


-V>br»  where  </>BR  satisfies 


k*  X  Av* 
sin*BR  "  2k  a  2X*  9  2V 


(3) 


wherein  V  is  the  sound  velocity. 

This  condition,  defining  the  Bragg  angle  is  similar  to  the  one  used  for 

X-ray  diffraction  from  the  periodic  arrangement  of  atoms  in  a  crystal.  As 

a  consequence,  the  incident  light  beam  is  deflected  from  its  original  direction 

through  the  angle  2^no.  If  we  introduce  angle  0  for  the  direction  of  propa- 

ok  n 

gation  of  the  n-th  order  light  beam,  then  from  (2)  we  obtain 


provided  all  angles  are  sufficiently  small. 

The  reader  should  be  warned  that  the  labelling  of  the  diffraction  orders 
throughout  this  paper  is  consistent  with  the  one  in  Raman-Nath's  original 
work  [42,43,48-52].  Hence,  the  positive  label  n  will  refer  to  the  diffracted 
light  beam  which  is  downshifted  in  frequency  by  an  amount  nuj*,  as  in  (l). 

In  theoretical  as  well  as  technical  papers,  written  by  (electrical)  engineers, 
the  positive  phasor  convention  is  used  [25],  consequently  the  positive  label 
n  would  refer  to  the  frequency  upshifted  light  beam  of  order  n  (ui^^w+nw* 
would  replace  (1)). 

Anyway,  the  frequencies  and  directions  of  propagation  of  the  diffracted  light 
waves  are  correctly  and  simply  determined  by  (1)  and  (4),  leaving  only  the 
amplitudes  and  thereafter  the  intensities  of  the  diffraction  orders  to  be 
calculated. 

In  a  series  of  extraordinary  papers  [42,43,48-52],  Raman  and  Nagendra  Nath 
use  at  first  geometrical  arguments  and  lateron  a  differential  equation  method 
to  determine  the  amplitudes  of  the  diffracted  light  beams.  In  their  elemen¬ 
tary  theory,  Raman  and  Nagendra  Nath  [48]  assumed  the  diffraction  effect  to 
be  due  to  a  simple  phase  corrugation  of  the  incident  light  caused  by  an 
acoustic  variation  of  the  index  of  refraction.  For  engineers,  phase  modulated 
waves  are  known  to  be  composed  of  a  carrier  and  a  set  of  sidebands  separated 
by  at*,  the  amplitudes  of  which  are  the  Besselfunctions  of  the  first  kind 
and  integral  order  n.  In  mathematical  terms,  the  modulated  output  electric 
field  can  be  represented  as 

E(x,z,t)  =  x(x,z,t)  exp  i(u)t-kz)  ,  (5) 

assuming  normal  incidence  for  the  input  light.  If  we  restrict  ourselves  to 
the  case  of  a  simple  sinusoidal  sound  field,  the  dielectric  permittivity  of 
the  disturbed  medium  can  be  written  as 

e(x,t)  ■  cgter  +  Ei  sin(w*t-k*x)]  .  (6) 

where  Cq  is  the  dielectric  constant  of  vacuum  and  is  the  peak  variation 
of  the  relative  permittivity  of  the  medium. 

For  pure  phase  modulation  due  to  the  progressive  sound  field  (6),  x(x,z,t) 
can  be  expressed  as  an  exponential  function  with  purely  imaginary  argument, 
i.e. 


r»\r*  i.  >  .  *  l 


ike .  z 

X(x,z,t)  =  exp[-^ —  sin(k*x-ui*t)]  . 

r 

Using  the  well-known  formula  of  Jacobi  [54] 


exp[(±i  cosx)z] 


(±i)nJ  (z)  exp(inx)  , 


which  generates  the  Bessel  functions  of  first  kind  and  integer  order  n, 
(7)  becomes 


n  ke .  z 

x(x,z,t)  *  £  Jn^2e — ^  exp[in(k*x-w*t)] 


Substituting  this  result  into  (5)  yields  the  Fourier  expansion  of  the 
electric  field  in  the  exit  plane  z=L  : 

+o» 

E(x,L,t)  **  ^  Jn(v)  exPCi(nkv'x-kL)]exp  i[(w-nu)*)t]  .  ( 10 ] 

n*“® 

In  (10),  the  constant 
ke.L  nc.L 

*-2r-  =  rf  <*»: 

r  r 

is  known  as  the  Raman-Nath  parameter  (RN  parameter)  and  it  is  a  measure  of 
the  strength  of  the  acoustic  modulation  of  the  light. 

From  (10)  we  learn  that  the  amplitude  of  the  n-th  order  diffracted  light 

beam  simply  is  J  (v);  hence  the  diffraction  order  irradiance  I  is 
n  n 


I  (v)  «  J„(v) 
n  n 


n  **  0,±1,±2, ...  . 


Raman  and  Nagendra  Nath  [49]  later  generalized  their  phase  modulation  model 
to  account  for  oblique  incidence.  Without  conceptual  changes  one  can  straight¬ 
forwardly  obtain  the  closed  form  expression  for  the  intensity  of  the  n-th 
order  diffracted  light  beam 

I  (v* )  ■  J^[v*  sinc(H(p^)3  ,  (13) 

n  n  AK 


with 


cos yp  2c  cos <fi  t  Xcosy 
r  r 


and 


,  VVVWVV'.'YTV  /V  V*  ^ 7JX  TUT  V-  ** »*  VT  v  unman  '  -  "  ~  ~ 


.  sin  nx 

sine  x  *  - 

TtX 


In  a  broader  perspective,  Raman  and  Nagendra  Nath  [42,51,52]  no  longer- 
restricted  the  diffraction  effect  to  be  exclusively  due  to  phase  modulation 
but  took  also  amplitude  modulation  into  account.  This  led  to  their  so-called 
generalized  theory.  Let  us  give  a  synopsis  of  this  since  then  well  developed 
approach. 

For  oblique  incidence,  instead  of  (5)  we  now  will  have 


E(x,z,t)  «*  x(x,z,t)  exp  i[ut-k(x  sin^»  +  z  cos»p)] 


where  the  amplitude  factor  x(x,z,t)  is  no  longer  restricted  to  the  form  (7), 
but  still  must  show  the  same  periodic  structure  as  the  disturbing  sound  wave. 
Consequently,  by  Fourier  expansion  the  function  x(x,z,t)  can  be  decomposed 
into  a  sum  of  plane  wave 


x(x,z,t)  ■  ^  $n(z)  exp  in(k*x*<u*t)  , 


hence. 


[x,z,t)  -  ^  $n(z)  exP  sinvOx  -  kz  cos?  ]  exp  i(u>-nu)*)t  .  (18) 


From  (18)  it  is  easy  to  rederive  (1)  and  (2),  whereas  the  unknown  amplitudes 


$n(z)  will  follow  from  inserting  the  expansion  into  the  appropriate  wave 


equation  for  the  medium  with  slowly  varying  permittivity  (6),  i.e. 


2  2  2 

3*E  .  3ZE  ,  3  l 

2  2  3  u0E(x.t)  -  • 

3x  3z  3t 


This  substitution  results  in  the  celebrated  Raman-Nath  equations  (called  from 
now  on  RN  equations), 


"  ^n-l"^n+l ^  “  in^nP’2a  sinv>)$n  ,  n  -  0,il,±2, 


linking  the  amplitudes  of  neighbouring  diffraction  orders  in  a  difference- 
differential  form.  Expressing  that  there  is  no  loss  of  energy  in  this  AO 
diffraction  effect,  provides  system  (20)  with  the  appropriate  (normalized) 
boundary  conditions 


♦  (0) 
n 


n  »  0,il ,i2,  . 


6 


In  (20),  we  introduced  the  phase  variable 


2t  ^  cosv>  L  8 


the  diffraction  regime  parameter 


*  ^  • 

the  latter  two  containing  the  ratio  of  wavelengths  of  light  and  ultrasound. 

For  a  physical  interpretation  of  the  parameter  a,  we  define  the  Bragg  angle 
of  higher  order  in  the  familiar  way  by 

Sin  *BR)  "  fir  “  2X^  "  •  n  “  0,±1,±2 .  (25) 

Then,  (20)  can  be  rewritten  as 

$  -  <WW  -  i"2»“  - •  <26> 

sin  *  BR 

Hence,  it  is  clear  that  a  sin^  is  a  measure  of  the  incidence  angle  in  units 
of  the  Bragg  angle 

After  Klein  and  Cook  [21],  quite  often,  the  parameter 


Q  *  pv  =*  ,  (27) 

A*Z 

is  introduced  to  devide  between  Bragg  and  Raman-Nath  diffraction  regimes. 

Since  this  parameter  is  independent  of  the  amplitude  of  the  disturbing 
sound  wave,  Q  is  necessary  but  not  sufficient  for  an  accurate  distinction 
between  both  regimes  [12,15,16,17,18,30,40,41,47,53].  Furthermore,  scrutini¬ 
zing  the  paper  of  Klein  and  Cook  reveals,  that  one  should  use  [15,16,17,29,53] 


i'  a  -2 — 


rather  than  Q  to  deal  with  cases  of  oblique  incidence. 


In  our  previous  work,  we  mainly  used  p  as  a  regime  parameter,  since  this 
one  gives  a  measure  of  the  amplitude  grating  impressed  on  the  light  wave. 
Indeed,  returning  to  (20),  for  p  =  0  we  obtain  the  well-known  Besselfunttion 
solution 


♦  n(C')  “  exp(-inbt')  , 


b  =  -  a  sim>  . 


Under  the  supplementary  restriction  of  normal  incidence  (^=0),  (26)  reduces 


*  (C)  =  J  (C)  , 

n  n 


kCjZ  vz  z 


From  (29),  resp.  (31),  at  z=L  we  rederive  the  results  (13),  resp.  (12),  as 
obtained  from  RN's  early  phase  modulation  model. 

For  normal  incidence,  we  obviously  have  C'  =  C,  V  =  v  and  Q'  =  Q. 

Let  us  remark  that,  recently,  there  was  a  revived  interest  in  equations  of 
the  RN  type,  not  only  for  their  relevance  to  acousto-optical  problems  [27] 
and  holography  [14,28],  but  also  to  the  study  of  other  physical  phenomena, 
such  as  stimulated  Compton  scattering  [6],  stimulated  Cherenkov  emission, 
and  to  the  theory  of  the  free  electron  laser  [9], 

In  the  following  we  shall  restrict  ourselves  to  normal  incidence  of  the  light, 
i.e.  ■  0,  so  that  the  system  of  RN  equations  reads 


n  2 

2dr*Un-r*n+l)"inp*n  *  n  "  °*±l*±2 . 


with  boundary  conditions 


♦„<°> 


n  ■  0,±1,±2,  ... 


It  has  been  shown  [32,38]  in  the  case  of  normal  incidence,  that 


w  v  www.  VM  to  a  ^■axoww?wm^> •>  "  '  ' 

8 


S  m-  W- 


so  that  the  system  (33)  can  be  simplified  into 
d*- 


df+*l  =  0  * 
d$  ? 

2dT  '  (Vl‘*n+l>  “  *"  H 


n  =  1,2,.. 


with  boundary  conditions 


(36) 

(37) 


*n(0)  =  6n0  *  n  "  O’1'2*  •••  •  (38) 

Eq.  (32)  has  as  immediate  consequence  that 

I_ n(0  =  In(0  ,  (39) 

meaning  that  for  normal  incidence  of  the  light,  the  diffraction  pattern  is 
symmetric  with  respect  to  the  zero  order. 

The  latter  system  has  been  solved  exactly  by  a  method  called  the  "generating 

function  method"  (GFM)  [26,37],  leading  to  an  infinite  series  expansion  for 

the  amplitudes  4>n(0  containing  the  Fourier  coefficients  of  the  even  periodic 

Mathieu  functions  ce2n(x,q).  Later  an  equivalent  method,  the  "modified 

generating  function  method"  (MGFM)  [15,16,17,34,35,36]  has  been  deviced, 

starting  directly  from  the  wave  equation  (19),  without  passing  through  the 

RN  equations,  and  giving  the  same  solution.  The  exact  solution  obtained  by 

those  methods  does  not  give  expressions  in  closed  form,  and  in  order  to 

perform  numerical  calculations,  one  should  have  the  disposal  of  extensive 

|  tables  of  the  Fourier  coefficients  of  the  periodic  Mathieu  functions.  The 

!  existing  tables  of  those  functions  are  either  too  restricted  [1]  or  are 

I  ^ 

related  to  Mathieu  functions  which  are  not  yet  expressed  in  the  canonical 

form  [7]  used  nowadays  [2,31].  This  is  the  reason  why  in  the  present  study 

we  shall  pay  close  attention  to  approximate  and  numerical  solutions  of  the 

system  (36,37)  with  boundary  conditions  (38). 

In  Section  1  we  give  a  survey  of  some  approximate  methods .  Section  2  considers 
in  detail  one  of  those  methods  (the  N-th  order  approximation  method)  which, 
reduced  to  an  eigenvalue  problem,  is  especially  suited  for  numerical  treatment 
In  the  next  Section  the  results  are  compared  with  experimental  data  of  Klein 
and  Hiedemann  [20,22]  and  in  Section  4  comparison  is  made  with  other  approxi¬ 
mations.  Finally  in  Section  5  some  further  considerations  with  respect  to 
the  N-th  order  approximation  method  are  expounded. 


.  -■>  .■>. 


> 

F, 
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NOMENCLATURE 


VU0 


permittivity,  magnetic  permeability  of  vacuum 


A,k,v  ,w,u 


wave  length,  wave  number,  frequency,  angular  frequency, 
speed  of  light  (in  medium) 


A*,k*,v*,u)*,V  wave  length,  wave  number,  frequency,  angular  frequency, 
speed  of  sound  (in  medium) 


e(x,t) 


J„U) 


MGFM 


dielectric  permittivity  of  the  disturbed  medium 


relative  permittivity  (undisturbed  medium) 


maximum  variation  of  the  linear  relative  permittivity 
(disturbed  medium) 


amplitude  of  the  (scattered)  electric  field  (disturbed  medium) 


acousto-optic  interaction  length  along  the  z-axis 


direction  of  sound  propagation 


direction  of  light  propagation  (at  normal  incidence) 


angle  of  incidence  of  light  (in  medium) 


order  of  diffraction 


Bragg  angle  of  order  n  (in  medium) 


deflection  angle  of  order  n  (in  medium) 


amplitude  (complex)  of  diffracted  light  wave  of  order  n 


complex  conjugate  of  <j> 


intensity  of  spectral  line  of  order  n 


Bessel  function  of  the  first  kind  of  order  n 


Raman-Nath  parameter  or  peak  phase  deviation 


regime  parameter 


Klein-Cook  parameter 


acousto-optic 


Raman-Nath 


generating  function  method 


modified  generating  function  method 


- .  v.\\,vv<v  ’>  ’> ' 


1.  APPROXIMATE  SOLUTIONS  OF  THE  RAMAN -NATH  EQUATIONS  (36.3?; 


1.1.  Raman-Nath's  elementary  theory 


The  most  simple  approximation  is  obtained  by  putting  p  **  0.  In  this  case  the 
infinite  system  (36,37)  with  boundary  conditions  (38),  or  the  equivalent 
system  (33)  with  boundary  conditions  (34),  has  an  exact  solution  (cf.  Eq.(31)) 


<J>n( v )  =  Jn(v)  ,  n  =  0,±1,±2,  ...  ,  (40) 

giving  the  intensities 

yv)  -  I.n(v)  *  J*(v)  ,  n  3  0,1,2,  ...  .  (41) 

A  proof  based  on  the  general  solution  of  the  system  (33)  with  p  3  0  is  given 
in  Appendix  A. 

A  comparison  of  this  elementary  theory  with  experiment  was  given  by  Klein 
[20]  and  by  Klein  and  Hiedemann  [22]  where  zeroth  and  first  order  light  inten¬ 
sities  vs  v  were  plotted  together  with  corresponding  intensity  measurements 
for  Q=0.94  (L=3.0cm),  1.26  (L=4.0cm)  and  1.48  (L=4.7cm)  (cf.  Figs.  5(p.), 

6(p. ) ,  7(p. ) ,  8(p. ) ,  9(p. ) ,  10(p. )). 

The  approximate  theory  gives  acceptable  results  only  if  v  *  2  for  Q=0.94, 
v  *  1.5  for  Q=1.26  and  v  =  1.3  for  Q=1.48.  Those  conditions  are  in  agreement 
with  the  theoretical  limits  established  by  Extermann  and  Wannier  [10]  as 
Q  «  2  and  Qv  <  2. 


1.2.  Perturbation  method 

In  order  to  obtain  an  improvement  of  RN's  elementary  theory  a  perturbation 
method  was  developed  by  Mertens  [32],  putting 

♦  (O  -  J  (C)  +  p*.(?)  +  pV,U)  +  ...  .  (42) 

n  n  nl  nz 

where  the  functions  ♦  ^(C)  and  ^^(C)  could  be  expressed  as  converging  power 
series  in  C-  Since  4^(0  is  of  the  form  C  o  ^(C).  Miller  and  Hiedemann  [39] 
wrote  this  solution  as 


Later,  using  the  generating  function  method,  Kuliasko  et  al  [26]  were  able 
to  express  the  solution  up  to  the  second -order  in  p  in  a  closed  form,  leading 
to  the  compact  formula  for  the  intensities 

In(v)  -  Jj(v)  -  y^Q  p2v2[(7n2-9v2)J2(v) 

\ 

+  (16n2-3)vJ  (v)J'(v)  -  5v2J,2(v)]  ,  n=0,± 1 .±2, . . .  ,  (44) 

n  n  n 

dJn(o 

where  J’(v)  -  -Jj—  >v  . 

One  of  the  principal  aims  of  Klein  [2Q]and  Klein  and  Hiedemann  [22]  was  the 
comparison  of  their  measured  light  intensities  with  the  intensities  obtained 
by  the  present  perturbation  theory.  For  Q=0.94  there  was  satisfactory  agreement 
between  experiment  and  theory  for  values  of  v  up  to  6.5  (cf.  Figs.  5(p.), 
8(p.));  for  Q=1.26  there  was  agreement  for  v  *  4.5  (cf.  Figs.  6(p.),  9 ( p . ) ) ; 
and  for  Q=1.48  the  agreement  was  rather  poor,  particularly  for  values  of  v 
larger  than  3.5  (cf.  Figs.  7(p.),  10 (p.)).  Klein  and  Hiedemann  conclude  that 
the  theory  gives  acceptable  results  for  experimental  arrangements  described 
by  Qv  <  8  and  that  an  upper  limit  of  usefulness  of  the  theory,  allowing  a 
reasonable  range  of  values  of  v,  is  Q  <  1.25.  Anyhow,  the  perturbation  method 
means  a  substantial  improvement  with  respect  to  the  results  obtained  by  the 
elementary  RN  theory. 


1.3. 


We  shall  now  discuss  a  totally  different  approximate  method  for  solving  the 
system  (36,37);  it  will  form  the  basis  of  our  numerical  treatment  of  the 
problem.  In  this  framework  we  neglect  the  energy  in  the  diffraction  orders 
higher  than  N,  i.e.  $+(N+1)  =  ^±(^+2)  “  ***  “  so  that  the  infinite  system 

(36,37)  is  replaced  by  the  following  truncated  system  of  N+l  equations  : 


with  the  boundary  conditions 


12 


$  (0)  «  6  n  ,  n  *  0,1,2, ...,N  . 
n  nu 


(46) 


This  procedure,  called  N-th  order  approximation  (NOA)  method  was  introduced 
by  Nagendra  Nath  [43]  for  N=1  (and,  hence,  restricted  to  the  Bragg  diffraction 
regime).  It  was  extended  for  arbitrary  N  by  Mertens  [33]  in  1962  for  the 
more  involved  problem  of  superposed  ultrasonic  waves.  In  that  paper  the 
denomination  NOA  method  was  coined.  Numerically,  however,  the  method  had  to 
be  restricted  to  N=2  (i.e.  five  equations),  according  to  the  computer  facilities 
at  that  time.  In  two  recent  publications  [4,5],  Blomme  and  Leroy  derived  finite 
analytical  expressions  for  the  intensities  using  respectively  N=2  and  N=3. 

A  profound  discussion  of  the  accuracy  of  those  20A  and  30A  methods,  in  a  wide 
range  of  the  parameters  p  and  v,  has  been  given  in  those  papers. 

To  end  this  Section  we  give  the  formulae  obtained  by  Nagendra  Nath  for  N=l, 
expressed  in  p  or  equivalently  in  Q  (=pv) 


J0(v)  "  1  ‘  “2 


2  v 


p  +8 


sin  ^  /p2+8  *»  1 


8v  2  /Q2+8v2 

2  sln  4  • 


Q2+8v 


2  v 


(47) 


4v 


1  (v)  »  __  sin  -  V+8  =  2 

11  pZ+8  Q  t8v 


sin 


2  /Q2+8v2 


A  derivation  of  those  formulae,  using  one  of  the  methods  exposed  in  the  next 
Section  is  given  in  Appendix  B. 

Those  expressions  will  be  usefyl  later  for  comparison  with  the  numerical 
results.  At  that  place  their  validity  shall  also  be  discussed. 


2.  SOLUTION  OF  THE  TRUNCATED  SYSTEM  (45)  WITH  BOUNDARY  CONDITIONS  (46) 

To  integrate  truncated  RN  systems  several  methods  have  been  proposed,  e.g. 
an  operator  technique  introduced  by  Benlarbi  and  Solymar  [3]  (extremely 
useful  to  treat  higher-order  Bragg  diffraction)  and  a  Laplace-transform 
method  [27].  In  this  paper  we  will  use  two  straightforward  techniques  : 

(1)  an  eigenvalue  method,  leading  to  a  characteristic  equation  of  degree  N+l, 
with  real  roots; 

(2)  the  operational  method  of  Heaviside-Jef freys  [19],  leading  to  expressions 
for  the  amplitudes  in  terms  of  determinants. 


Assuming  a  solution  of  system  (45)  in  the  natural  form 
♦  »  an  exp(^isC) 

the  corresponding  eigenvalue  problem  would  have  a  non-Hermitian  matrix,  the 
condition  for  real  eigenvalues  thus  being  abolished.  This  slight  problem, 
which  is  due  to  the  use  of  symmetry  relation  (35),  can  be  easily  overcome 
by  introducing  a  similarity  transformation,  realized  by  directly  modifying 
the  solution  into 

♦0  *  ^2  aQ  exp(|is; )  , 

(48) 

♦n  =  aR  exp(^-isC)  •  n  =  1,2,...,N  . 

Substitution  of  (48)  into  system  (45)  leads  to  a  system  of  linear  homogeneous 
equations  for  a^,  ,a^ , . . .  ,a^,  which  may  be  written  in  matrix  form  as 

(M-sI).a  -  O  ,  (49) 


T 

where  X  is  the  (N+1)*(N+1)  unit  matrix,  a  =  (a_  a,  ...  a.,)  and 

u  l  N 


is  a  (N+1)*(N+1)  Hermitian  matrix,  the  eigenvalues  of  which  are  real  [11]. 

In  order  to  obtain  a  nonzero  solution  a  of  the  linear  homogeneous  system  (49), 


the  determinant  of  the  coefficients  must  vanish  : 


det(M-sI)  =  0  . 


In  other  words,  we  ought  to  calculate  the  N+l  real  eigenvalues  s.  (k-1,2,..., 

/  I  \  K 

N+l)  of  tn^  Hermitian  matrix  M.  The  eigenvector  a.'  ,  with  a.(k)T  * 

(k)  (k)  (k) 

(a«  a,  ...  a„  ),  associated  with  the  eigenvalue  s.  will  then  be  determined 
U  1  N  K 

by  the  linear  homogeneous  system 

(M-skI).a(k)  =  O  .  (52) 


(k)  (k) 

From  the  structure  of  the  system  (52),  we  may  see  that  a^  ^  0  and  that  a^ 

may  be  chosen  real  if  n  is  even,  and  purely  imaginary  if  n  is  odd. 

(k)  — 

Furthermore,  we  can  choose  a^  «  -Jlil  (k=»l,2 . N+l)  since  any  eigenvector 

is  only  determined  up  to  an  arbitrary  factor.  The  general  solution  of  the 
truncated  linear  differential  system  (45)  may  then  be  written  as  [11] 

N+l 

♦o  *  I  ck  e’‘P(2iskU  •  (53) 

k=l 

N+l 

♦n  *■  £  ckank^  exP^2iskC’^  ’  n  =  1,2,...,N  ,  (54) 

k=l 

(k) 

uhprp  t"ViP  N+1  rnnct- anf «:  r  are  rpal  ( a<s  a  rnnspmipncp  of  ths  choices  of  a_ 


where  the  N+l  constants 

c,. 

and  a^k^).  These 
n 

JV 

constants 

to  (53)  and  (54) 

• 

• 

N+l 

r* 

N+l 

L  cx*‘ 

V 

’  L 

ck 

k=l 

k=l 

Finally  we  can  calculate  the  intensities  in  z«L  : 

N+l 

Iq(v)  »  0q(v)$o(v)  =  1  -  4  £  c^ck  sin2(Sj-sk)^  , 

j.k-1 

j<k 


I.„(v)  "  “  "4  )  c.c.a^V 

in  in  in  /t  j  k  n  n 


(k)  .  2,  xV 

„  sin 


j,k=l 

j<k 

where  the  bar  stands  for  complex  conjugate. 


n  ■  1 , . .  . ,  M  , 


VVi 


Needless  to  say  that  the  characteristic  equation  (51),  of  degree  N+l  in  s, 
can  only  be  solved  analytically  for  N43.  Thus  only  in  the  latter  cases 
explicit  analytical  expressions  for  the  intensities  can  be  obtained  [4,5,33,43] 
Otherwise  the  problem  has  to  be  treated  numerically  using  the  following  steps  : 

(i)  determine  the  eigenvalues  and  eigenvectors  of  the  matrix  M; 

(ii)  next  solve  the  linear  system  (55)  for  c^; 

(iii)  substitute  the  results  obtained  from  (i)  and  (ii)  into  Eqs.  (56)  and 
(57). 

We  remark  that  the  solution  scheme  of  the  eigenvalue  problem  sketched  above 
is  particularly  suited  to  computer  algorithms. 


2.2.  Heaviside's  operational  method 


Now  we  will  apply  Heaviside's  operational  method  to  the  truncated  system  (45). 
After  Jeffreys  [19,  p.237]  we  write  p  for  d/dC  and  interpret  p  ^  as  the 
operation  of  definite  integration 

p  1f(C)  =  [  f (z)dz  .  (58) 

J0 

The  resulting  subsidiary  equations  [19],  which  take  the  boundary  conditions 
(46)  into  account,  are 


-s0Q  +  2i01  =  -s 

(pn2-s)<jin  =  0  ,  n  «  1,2 . N  ,  (59) 

(pN2-s)$N  -  id»N-1  *  0 

with  s  ■  -2ip.  This  system  is  now  to  be  solved  as  if  p  (or  s)  were  a  number. 
In  compact  notation  the  system  (59)  can  be  written  as 


D<t>  ■  -sE  , 


(60) 


where  the  almost  Hermitian  (N+l ) * (N+l )  matrix  D  is  related  to  M  through 
a  similarity  transformation 

D(s)  -  P(M-sX)P'1  ,  (61) 

with  P  ■  I+R,  such  that  all  elements  r^  of  R  vanish,  except  r^  *»  /2-1. 


T  T 

Furthermore  we  have  <t>  =  ($q  $2  •••  and  E  =  (1  0  ...  0),  the  latter 

vector  expressing  the  boundary  condition.  From  (60)  we  obtain  the  formal 
solution  of  the  problem  : 

<t>  -  -sD_1E  .  (62) 

where  ID  *  is  the  inverse  of  3D.  So,  one  explicitly  has 


1  ,n+l 
det  ID 


n  0,1, ...,N  , 


where  D.  stands  for  the  cofactor  of  the  element  d,  . ,  of  3D  (n=0,l,2, . . . ,N) 
l,n+l  l,n+l 

Substituting  s  =  -2ip  in  (63)  we  find 


2pFl,n+l(p)  ,  ,  2p 

*n  '  “Hjj -  (n=1M . N'U  •  *N  '  F(57  • 


where  Fj  n+^(p)  is  the  determinant 


2p-i(n+l )  p 


2p-i(n+2)  p  1  0  ... 


0  -1  2p-i(N-l)  p 


0 


2p-iN2p 


and  F(p)  ■  (-1)”  aet  ID(p). 

Next,  in  order  to  express  formulae  (64)  in  the  original  variable  C,  we  will 

apply  Heaviside's  expansion  theorem  [19,  p.238],  i.e. 

N+l  . 

A(p)  _  A(0)  .  V  MV 


A(p)  A(0)  .  V  AVV  , 

S5T  mm  +  l  a/B’C.)  • 


where  A(p)  is  a  polynomial  in  p  of  the  same  or  a  lower  degree  than  B(p); 
0^  (k«l , . . . ,N+1 )  are  the  simple  zeros  of  B(p)  and  B'  ■  dB/dp.  Hence,  (64) 
can  be  replaced  by 
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N?‘  Fi  „+,<V 

♦  n  -  2  |  — pi  —  exp(ak^)  ,  n  =>  0,1,...,N-1  ; 

k-1  k 


exp(ak;) 


F*(ak)  * 


Since  from  (61),  the  matrices  D  and  M-sI  are  similar,  one  can  prove  that 


ok  ■  2^Sk’  w^ere  sk  rePresent  the  real  eigenvalues  of  M.  After  some  lengthy 


calculations  we  obtain  the  expressions  for  the  intensities  of  the  spectral 


lines 


r  F  (yis.)F..( 
I„(v)  -1-16  l  ljk 

j  ,k=l  F'(jis^.)F'(^isk) 


Fll(2isi)Fll(21Sk)  .  2,  _  vv 

-  sin  \S.“S,  ),  , 

A  .  .  A  .  JKA 


[  (v)  -  -16  )  -AfUtL?— 3 — i 

n  L  l  * 

j,k*l  F ' (yis . )F 1 (ris 

1  /l.  J 


(yis.)F.  .(yis,)  - 

2  i  l,n+l  2k  .2,  \V 

—=7- -  Sln  Sj*Sk  4 


n  «  1,2,...,N  ,  with  F,  ..,.(ks.)  =  1,  and  where  we  have  taken  into  account 

1 »N+1  2  j  . 


that  the  determinants  F ,  ,  .(lis.)  and  F(ris.)  are  either  real  or  purely 

1  ,n+l  2  J  2  j  v 

imaginary.  In  order  to  simplify  the  calculations  of  the  determinants  it  may 
be  advantageous  to  use  the  following  recursion  relations  ; 


VP>  -  <2p-ij  P)F1>J+1(p)  +  Fl  j+2(p)  ,  Fj_n+1  -  1  ,  F 


1  ,N+2 


F'(p)  -  2Fu(p)  +  2pF|1(p)  +  2Fj2(p)  s 


Fij(p)  *  2Fl,j+l(p)  +  (2p-ij2p)F|t j+1(p)  +  F|>  ,+2(p)  ,  j  -  1 . N 


It  is  clear,  especially  regarding  (61),  that  the  eigenvalue  method  (Subsection 
2.1)  and  Heaviside's  operational  method  (Subsection  2.2)  are  completely 
equivalent. 

At  this  point  it  is  perhaps  worthwhile  to  comment  on  the  paper  by  Leroy  and 

Claeys  [27]  wherein  a  Laplace  transform  method  is  used  to  integrate  the 

truncated  system  (45).  To  begin  with  the  Laplace  transform  F(s)  •  £[f(z)]  * 

•  ^2 

/o  f(z)e  dz  imposes  the  function  f(z)  to  be  defined  in  a  domain  [0,<»[  for 
the  variable  z.  However,  the  functions  $n({;)i  with  C  "  uc^z/Aer,  describe 
the  amplitudes  of  the  diffracted  light  in  the  finite  domain  [0,v],  where  the 
ultrasonic  field  is  active.  So  without  either  changing  the  variables,  extending 
the  domain  of  the  function  towards  [0,«[  or  introducing  a  Heaviside  unit 


f 


wv! 


function  restricting  the  domain  to  [Q,L],  the  use  of  a  Laplace  transform  in 
the  present  form  is  physically  inadmissible.  To  support  this  case  we  refer 
to  Jeffrey's  comments  [19,  pp. 458, 459]  on  the  Laplace  transform,  where  next 
to  physical  arguments,  mathematical  plea  in  favour  of  Heaviside's  operational 
method  is  given.  Secondly  the  calculation  of  the  inverse  transforms  $  (^) 
in  Leroy  and  Claeys'  paper  leads  to  rather  complicated  infinite  series  of 
Bessel  functions,  contrary  to  the  finite  expressions  for  the  amplitudes  obtained 
by  Heaviside's  expansion  theorem  used  above. 


3.  COMPARISON  OF  THE  DIFFERENT  APPROXIMATE  METHODS  WITH  EXPERIMENTAL  RESULTS 

In  Fig.  5  (Q=0.94),  Fig.  6  (Q=1.26)  and  Fig.  7  (Q=1.48)  we  compare  curves 
for  Iq  vs  v  obtained  from  RN's  geometrical  theory  [48],  Mertens*  perturbation 
theory  [26,32]  and  the  NOA  method  (for  N=7 )  with  the  experimental  results 
of  Klein  and  Hiedemann  [22].  In  Fig.  8  (Q=0.94),  Fig.  9  (Q*1.26)  and  Fig.  10 
(Q=1.48)  curves  for  Ij  vs  v  obtained  with  the  same  approximate  methods  are 
compared  with  the  experimental  points  obtained  by  Klein  [20].  It  must  be 
remarked  that  in  the  experiments  of  Klein,  according  to  a  finite-amplitude 
distortion  of  the  sound  wave  an  asymmetric  diffraction  pattern  was  observed. 

In  all  plots  of  1^  the  average  of  the  measured  positive  and  negative  light 
intensities  is  shown  and  in  no  case  the  measured  intensities  differ  at  most 
5  percent.  From  numerical  calculation,  it  was  found  that,  for  the  considered 
values  of  Q,  a  0  for  j  2  7;  a  typical  example  is  shown  in  Fig.  11  for 
Q-0.94.  Hence,  we  have  computed  the  intensities  (56),  (57)  and  (68),  (69) 
up  to  N=*7.  The  theoretical  results  perfectly  coincide  since  both  methods 
considered  are  equivalent;  not  surprisingly  concerning  CPU  time,  calculations 
with  Heaviside's  operational  method  are  about  25%  faster  than  the  ones  based 
on  the  eigenvalue  method.  For  the  three  considered  values  of  Q  the  theore¬ 
tical  curves  fit  the  experimental  points  for  1^  and  1^  perfectly. 

As  discussed  earlier  (Subsection  1.1)  it  is  clearly  seen  from  the  figures 
that  the  deviation  from  RN's  elementary  theory  is  already  apparant  for  rather 
low  values  of  v.  Evidently  there  is  an  excellent  fitting  with  the  experimental 
points  whereas  the  perturbation  method  fails  (cf.  the  discussion  in 
Subsection  1.2). 
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4.  COMPARISON  OF  THE  DIFFERENT  APPROXIMATE  METHODS  TO  ONE  ANOTHER 


In  this  section  the  NOA  method,  with  N  ranging  from  7  to  15,  is  carried  out 
for  selected  values  of  the  Klein-Cook  parameter  Q  (between  0.1  and  50)  and 
with  the  RN  parameter  v  ranging  from  0  to  15.  The  zeroth  and  first  order 
intensities  thus  obtained  are  compared  with  either  the  squared  Bessel  function 
expressions  or  with  the  squared  sine  expressions  obtained  from  the  10A  method 
(Eqs.  (47)). 

In  Fig.  12  we  compare  the  curves  for  I»(v)  and  l.(v),  computed  from  (56)  and 

(57)  (n=l)  for  Q=0.1  and  Q=l,  with  the  squared  Bessel  functions  J"(v)  sn<* 

2  U 
J^(v).  Remark  that  for  Q=0.1,  the  zeroth  and  first  order  intensities  perfectly 

coincide  with  the  corresponding  Bessel  function  expressions.  For  Q=5  and  Q=50, 

Iq(v)  and  I^(v)  are  respectively  shown  in  Figs.  13  and  14.  For  Q=0.1  it  was 

necessary  to  take  N=15,  i.e.  a  system  of  16  coupled  equations  had  to  be  taken 

into  account.  For  all  other  values  of  Q  a  value  of  N=7  was  sufficient.  To 

test  our  numerical  procedure  based  on  the  eigenvalue  method  we  have  compared 

our  results  with  similar  ones  computed  by  straightforward  numerical  integration 

of  the  system  (45)  by  a  discrete  step  method  [21,24,45],  From  Figs.  11  to  13 

(and  additional  figures  for  stepwise  values  of  Q  between  0.01  and  100,  not 

shown  here)  we  may  conclude  that 

-  For  Q<<1  (i.e.  Q  <  0.1)  the  intensities  Iq  and  1^  calculated  with  the  NOA 
method  coincide  practically  with  squared  Bessel  functions  for  the  whole 
range  of  v  considered; 

-  For  0.1  <  Q  <  2  the  intensities  fit  the  squared  Bessel  function  expressions 
up  to  a  threshold  in  the  close  neighbourhood  of  v=2.405  (the  first  zero 

of  Jq( v) ) ; 

-  For  2  <  Q  <  10  the  intensities  considerably  deviate  as  well  from  the 
squared  Bessel  functions  as  from  the  squared  sine  expressions  (for  Q=5 
there  is  only  a  good  fitting  with  the  latter  up  to  v  a  2); 

-  For  Q>>1  (i.e.  Q>10),  up  to  a  certain  threshold  v^  (v^  s  5  for  Q=50), 
practically  all  the  energy  of  the  incident  light  is  concentrated  in  Iq 
so  that  no  diffraction  occurs;  above  that  threshold  diffraction  only 
excites  first  orders  (Bragg  diffraction);  the  fitting  of  the  NOA  curves 
(N»7)  with  the  squared  sine  curves  (N=l)  is  excellent  up  to  2v^(  a  10  for 
Q-50). 


5.  FURTHER  NEW  RESULTS  OBTAINED  WITH  THE  NOA  METHOD 


As  argued  in  recent  work  [12,14,53],  p,  v  and  Q=pv  are  relevant  (though 
dependent)  parameters  for  the  investigation  of  criteria  for  multi-order  RN 
diffraction  or  two-order  Bragg  diffraction.  Based  on  our  numerical  integration 
method,  outlined  above,  we  present  in  Fig.  15  a  three-dimensional  plot  of 
Iq  as  function  of  Q  and  v  (similarly  one  could  use  p  and  v).  The  cross-sections 
correspond  to  constant  values  of  Q  ranging  from  0  to  10  with  steps  of  0.25. 
Analogously  in  Fig.  16  1^  is  plotted  as  a  function  of  Q  and  v.  In  both  figures 
one  may  clearly  observe  the  evolution  from  the  RN  regime  (Q<<1)  to  the  Bragg 
regime  (Q>>1),  passing  trough  all  intermediate  steps.  In  Fig.  17  we  show 
curves  for  constant  Iq  in  a  v,Q-plane.  In  order  to  realize  a  chosen  constant 
value  of  Iq,  this  plot  allows  the  selection  of  corresponding  characteristic 
couples  (v,Q).  Fig.  18  presents  an  analogous  plot  for  1..  In  Figs.  19,  20, 
curves  of  constant  Iq,  respectively  1^,  are  depicted  in  a  Q,p-plane,  from 
which  pairs  (Q,p)  may  be  obtained.  It  is  apparant  in  both  figures,  that  the 
behaviour  of  the  iso-curves  drastically  changes  when  passing  from  the  region 
where  v<4  to  the  one  where  v>4.  In  the  region  above  the  line  p  =  Q/4  (v<4) 
the  iso-curves  present  a  rather  simple  structure,  while  in  the  region  below 
the  line  p  =  Q/4  (v>4)  they  are  behaving  rather  chaotically.  A  more  profound 
analysis  and  interpretation  of  the  results  will  be  presented  later. 


APPENDIX  B 


An  example  of  the  use  of  the  eigenvalue  problem  developed  in  subsection  2.1 


We  shall  take  N=l,  which  will  lead  to  an  analytical  solution  of  the  problem. 
The  truncated  system  (45)  becomes 


-  +  A  =  0 

dC  U  ' 


2JT  '  *0  ’  ioh  • 


with  boundary  conditions 


$q(0)  =  1  ,  <^(0)  «  0 


As  a  solution  we  assume 


«>0  -  /2  aQ  exp(^isC)  , 


4^  -  exp(|isO  . 


Upon  substitution  in  Eqs.  (Bl)  the  characteristic  equation  follows 


det(M-sI)  = 


s  -  ps  -  2  =  0  , 


-i/2  p-s 


where  the  matrix 


-i/2  p 


is  Hermitian.  From  (B5)  the  (real)  eigenvalues  are  easily  obtained  : 


S1  »  2^p  +  ^p2+8)  » 


s2  -  2<p  -  / P:+8)  . 
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The  equation  for  the  eigenvalues  is 


(M-sRX).a.(k)  =  O  . 


For  k=l  we  thus  obtain 


+  i/2aJU  =  0  , 
-i/Ia^  +  (p-s.la^  =  0  , 


(BIO) 


where  the  second  equation  is  in  fact  redundant.  A  solution  may  be  chosen  as 


(Bll) 


1 

2  S1  * 


Similarly,  we  find  fod  k=2 


(2)  _  /2 


(B12) 


(2)  i 
ll  "  "  2  S2  * 


Hence,  the  amplitudes  can  be  written  as  (cf.  Eqs.  (53)  and  (54)) 


<frQ  ■  c1  exptjis^)  +  c2  exp(|is2C)  , 


(B13) 


*1  ■  c  exp(jis1C)  +  c2  a|2^  exp(|is2(;)  . 


(B14) 


The  constants  are  determined  by  the  boundary  conditions  (B2),  hence. 


C1  +  c2  "  1  • 


(B15) 


c  a(1)  +  c  a(2)  -  0 
clal  +  c2al  U 


Taking  into  account  (Bll)  and  (B12)  we  obtain 


c  -  -  - 

1  V5z 


•2  s,-s2  ' 


(  B 1 6 ) 


Calculating  the  intensities  from  the  amplitudes  (B13)  and  (B14),  therefore 
using  the  values  of  the  constants  and  c^.  we  find 


!0  ‘  Vo  "  1  +  ;  ‘  A;  si"2  4(VS2K  ' 


(B17) 


(srszr 


1tl  "  *±l*±l 


sls2  2  21 

("'  .")  sin^ 


S1'S2 


(B18) 


Finally,  introducing  the  explicit  expressions  for  the  eigenvalues  (Eqs.  (B7) 
and  (B8))  we  recover  Nagendra  Nath's  classical  results 


IQ  -  1  -  -|—  sin2  |v^+8  C  ,  (B19) 

p  +8- 

I,.  -  I  ,  -  -4—  sin2  rSP+8  C  • 


(B20) 
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SAMENVATTING 


De  bedoeling  van  onderhavige  studie  is  een  overzicht  te  geven  van  klassieke 
en  meer  recente  benaderde  methoden  in  het  akoesto-optisch  dif fractievraagstuk 
en  hun  numerieke  resultaten  te  vergelijken  met  de  experimented  gegevens  van 
Klein  en  Hiedemann. 

In  de  inleiding  wordt  de  diffractie  van  licht  door  ultrageluiden  beknopt 
beschreven  en  wordt  aangetoond,  dat  de  golfvergeli jking  van  het  licht,  dat 
zich  voortplant  in  een  vloeibaar  medium  door  een  ultrageluidsgolf  gestoord, 
aanleiding  geeft  tot  een  oneindig  stelsel  dif fero-dif f erentiaalvergelijkingen 
voor  de  amplituden  van  de  gedif fracteerde  lichtgolven;  dit  stelsel  is  het  zg. 
Raman-Nath  (RN)  stelsel,  waarvan  de  n  vergelijking  (n=0,±l ,±2, . . . )  het  verband 
geeft  tussen  de  amplituden  van  de  orden  n-1,  n  en  n+1;  deze  amplituden  voldoen 
aan  eenvoudige  randvoorwaarden.  Drie  parameters  spelen  hierbij  een  belang- 
rijke  rol 


(1)  de  zg.  Raman-Nath  parameter  v  =  nc.L/e  X  ;  (2)  de  regime- 

2  2  1  r  2 

parameter  p  =*  2e  \  /t^\*  ;  (3)  de  zg.  Klein-Cook  parameter  Q  =  pv  «*  2tiAL/A* 


Terwijl  de  eerste  twee  afhankelijk  zijn  van  e^,  die  een  maat  is  voor  de 
sterkte  van  de  ultrageluidsgolf,  is  de  derde  parameter  daarvan  onafhankelijk. 

In  een  eerste  paragraaf  wordt  een  overzicht  gegeven  van  drie  benaderde  methoden 
voor  de  behandeling  van  het  Raman-Nath  stelsel. 

(1)  p  *  0,  hetgeen  benaderd  verwezenlijkt  wordt  voor  grote  golflengten  van  de 
ultrageluiden.  In  dit  geval  kan  het  oneindig  RN  stelsel  exakt  worden  opgelost, 
hetgeen  leidt  tot  dezelfde  resultaten  als  deze  bekomen  met  de  elementaire 
geometrisch  optische  theorie  van  Raman  en  Nagendra  Nath.  De  intensiteiten  van 


de  verschillende  diffractielijnen  worden  in  dit  geval  gegeven  door  de 


eenvoudige  uitdrukkingen  I 


J  (v),  waarbij  J  (v)  de  Besselfuntie  van  orde  n 
n  J  n 


voorstelt.  Deze  benaderde  theorie  wordt  door  Klein  en  Hiedemann  met  hun 
experimentele  resultaten  vergeleken.  Ze  blijkt  aanvaardbaar  voor  v  <  2  als 
Q»0.94,  v  <  1.5  als  Q=1.26  en  v  <  1.3  als  Q=1.48. 

(2)  Een  verbetering  van  deze  elementaire  RN  theorie  wordt  bekomen  door  de 
storingstheorie,  die  vooreerst  door  Mertens  werd  opgesteld,  waarbij  de  ampli¬ 
tude  van  iedere  gedif fracteerde  golf  als  een  reeksontwikkeling  naar  p  (of  naar 


Q)  wordt  voorgesteld,  en  waarvan  de  eerste  term  Jn(v)  Is*  Vergeleken  met  de 


experimenten  van  Klein  en  Hiedemann  betekent  ze  een  verbetering  t.o.v.  de 
elementaire  RN  theorie.  De  overeenstemming  tussen  theorie  en  experiment  is 
echter  niet  bevredigend  voor  v  >  6.5  als  Q=0.94,  v  >  4.5  als  5=1.26  en  v  >  3.5 
als  Q-1.48. 


(3)  De  NOA  methode  (Ne  orde  approximate) .  Deze  methode  bestaat  er  in  het 


6 

RN  stelsel  te  benaderen  door  de  gediff racteerde  orden  hoger  dan  de  N  te 
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verwaarlozen,  zodat  het  gesymmetriseerde  RN  stelsel  nog  maar  N+l  vergelijkingen 
bevat.  Deze  methode  werd  voor  het  eerst  in  1939  door  Nagendra  Nath  voor  N=1 
opgesteld  en  later  door  Hertens  tot  willekeurige  waarden  van  N  uitgebreid. 

De  tweede  paragraaf  behandelt  de  oplossing  van  dit  eindige  (afgebrokeri)  stelsel 
differentiaalvergelijkingen  met  gegeven  randvoorwaarden.  Twee  methoden  worden 
gebezlgd  :  (1)  een  eigenwaardenmethode  die  leidt  tot  een  karakteristieke  verge- 
lijking  van  de  graad  N+l  met  reele  wortels  (eigenwaarder^ ;  de  intensiteiten 
worden  dan  uitgedrukt  als  een  eindige  son  van  kwadraten  van  sinussen,  die 
functie  zijn  van  v  en  van  de  verschillen  van  de  eigenwaarden  twee  aan  twee; 
de  coef f icienten  zijn  afhankelijk  van  de  componenten  van  de  eigenvectoren; 

(2)  de  operationele  methode  van  Heaviside- Jeff reys ,  die  de  amplituden  van  de 
gedif f racteerde  golven  uitdrukt  in  functie  van  determinanten. 

Beide  methoden,  die  trouwens  geli jkwaardig  zijn,  zijn  uitstekend  geschikt  voor 
de  behandeling  van  het  akoesto-optisch  vraagstuk  met  behulp  van  een  computer. 

De  vergelijking  van  de  NOA  methode  met  de  experimentele  gegevens  van  Klein 
en  Hiedemann  maakt  het  onderwerp  uit  van  de  derde  paragraaf.  Uit  de  computer- 
berekeningen  blijkt  dat  voor  de  beschouwde  waarden  van  Q  (0.94,  1.26,  1.48), 

N=7  een  uitstekende  benadering  geeft,  die  bovendien  buitengewoon  goed  bij  de 
experimentele  punten  aansluit. 

In  paragraaf  4  worden  de  verschillende  behandelde  benaderingsmethoden  syste- 
matisch  met  elkaar  vergeleken  voor  Q=0.1,  1,  5  en  50. 

In  een  laatste  paragraaf  tenslotte  worden  nog  enkele  nieuwe  resultaten  aange- 
haald,  bekomen  door  middel  van  de  NOA  methode,  nl., 

-  de  oppervlakken  van  Iq  en  1^  als  een  functie  van  Q  en  v; 

-  de  isokrommen  voor  Iq  en  1^  in  een  v,Q-diagram; 

-  de  isokrommen  voor  IQ  en  1^  in  een  Q,p-diagram. 

Grondige  interpretatie  van  deze  resultaten  en  het  eventueel  verband  met  een 
keuze  van  p  of  Q  voor  het  bepalen  van  het  dif f ractieregime  zal  nader  worden 
onderzocht . 
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